WHAT DOES QUANTUM CHAOS MEAN? 
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We discuss how the concept of the quantum action can be used to characterize quan- 
tum chaos. As an example we study quantum mechanics of the inverse square po- 
tential in order to test some questions related to quantum action. Quantum chaos 
is discussed for the 2-D hamiltonian system of harmonic oscillators with anharmonic 
coupling. 



The classical chaos has no direct analogue in quantum physics. One reason is that 
in quantum mechanics (Q.M.) one cannot define a point in phase space (due to 
Heisenberg's uncertainty relation), hence one cannot define Lyapunov exponents from 
diverging trajectories, thus quantum chaos cannot be represented quantitatively via 
Poincare sections. In Ref. (Cammetti et al., 2000) the effective action has been used to 
describe classically chaotic quantum systems. In Ref. (Jirari et al., 2001) the quantum 
action has been introduced. It basically states that transition amplitudes in quantum 
mechanics can be expressed in terms of an action, which has the form of the classical 
action, but with parameters (mass, potential) which are different from those of the 
classical action. The quantum action can be interpreted as a renormalized action in 
Q.M. (Jirari et al., 2002). It is similar to the standard effective action, however, it 
is free of the deseases of the latter (infinite series of higher derivative terms, non- 
localities etc.). In the limit of large imaginary transition time, the existence of the 
quantum action has been proven (Kroger, 2002). 

The quantum action has been originally proposed as a conjecture obeying follow- 
ing requirements. For a given classical action 
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1. INTRODUCTION 





there is a quantum action 




(2) 
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which parametrizes the quantum mechanical transition amplitude (Green's function) 
G(y, t, x, 0) by the expression 

G(y, t, x, 0) = Zexp -t(x d )/h , (3) 

where x c i denotes the classical path corresponding to the action S and 
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E(x cZ ) = S[x 
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where t is the transition time, x and y are final and initial points of trajectory and Z is 
normalization factor. A priori, it is not evident that such quantum action exists. By 
now the existence of the quantum action has been established in the following cases: 
(i) Harmonic oscillator. In this case the quantum action is identical to the classical 
action, (ii) In the limit when the transition time t — > 0, the quantum action exists. 
Dirac has noticed that the path integral in this limit is dominated by the contribution 
from the classical trajectory. Hence the quantum action coincides with the classical 
action, (iii) In imaginary time (t — > —iT is necessary to describe thermodynamics 
at finite temperature) and going to the limit T — > oo, the quantum action exists. In 
general it is quite different from the classical action. 

This leads to the question: What about finite transition time? Does the quantum 
action exist for arbitray finite transition time? Does it parametrize Q.M. transition 
amplitudes for all x to y equally well? In general, the quantum action has to be 
determined non-perturbatively by numerical computations. What can be said about 
stability, convergence and errors of such procedure? What is quantum chaos in the 
framework of the quantum action? Those questions are the subject of the paper. 

2. INVERSE SQUARE POTENTIAL 

We consider the following classical potential in 1-D 

V(x) = -mu 2 x 2 + g x~ 2 = v 2 x 2 + v -2 x~ 2 . (4) 

The corresponding classical action is given by S = J dt ^mx 2 — V(x), while the 
classical action in imaginary time (t — > —iT) is given by Se = J dt ^mx 2 + V(x) 
The potential is parity symmetric. Because it has an infinite barrier (for g > 0) at 
the origin, the system at x < is separated from the system at x > 0. We consider 
only the motion in the domain x > 0. The potential is shown in Fig. pp. We have 
chosen to consider the inverse square potential, because of the distinct feature that 
the corresponding quantum mechanical transition amplitudes are known analytically 
(Khandehar & Lavande, 1975; Schulman, 1981, sec. 32). 
The Euclidean transition amplitude reads 
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1 + 8mg/h ) / 2. The transition 
amplitude contains all information on the spectrum and wave functions. For example, 
by going to the limit T —>■ oo (Feynman-Kac limit) one reads off the ground state 
energy and wave function 

E gr = Ml + 7), lU*) = Zl 12 x 1 ' 2 ^ exp[-^x 2 ] . (6) 

This is, of course, identical with the direct solution from the Schrodinger equation. 
The wave function is shown in Fig. £Q . 

Wave function with classical and quantum 
potential for g=l/100 

\ Particle in the potential : V(x) = mco 2 x 2 /2+g/x 2 




Figure 1: Inverse square potential. Classical potential with (full line), quantum 
potential (dashed line) and ground state wave function ( dotted line). 

2.1. Dynamical time and length scales 

In this work we do a numerical study of the inverse square potentialal to test the 
quantum action. The numerical studies presented below have been done in imagi- 
nary time. Our choice of the model has been influenced by the analytical solvability 
and less by the question if it plays a role in nature. Consequently, absolute values 
in physical units of the model parameters (mass, potential parameters) are not of 
primary importance. We have expressed those parameters in dimensionless units. 
However, what is important are time and length scales, which are dynamically gen- 
erated by the model, say a time scale T sc and a length scale A sc . Those scales serve 
as reference values to give a sense to statements like "for large transition times T", 
or a "small spatial resolution Ax", which means T/T sc » 1 and Ax/A sc « 1, 



respectively. How to choose those dynamical scales? We have used as time scale 

Tg C \j Eg r . 

As a length scale one may introduce the analogue of the Bohr radius of the ground 
state wave function. Another possibility is to define a length scale A sc by 
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dx \^ gr {x)\ 2 = 0.95 , (7) 



i.e. the length which covers 95% of the probability of the ground state wave function. 
For the classical action considered in this paper those scale parameters are T sc = 0.4 
and A sc « 2.35. 

2.2. Analytical and numerical results for the quantum action in the asymp- 
totic regime 

In Ref.(Jirari et al., 2002) we have shown for the Euclidean asymptotic regime T — > oo, 
that the following analytic relations exist between ground state energy E gr , ground 
state wave function ip gr and the quantum action. V m in = E gr , i.e. the minimum of 
the quantum potential gives the ground state energy. Next, the ground state wave 
function can be expressed in terms of the parameters of the quantum action 



^gr{x) = — e 



1 o'Sl ■ dx ' \flrh{V{x')-V min )/h 
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Finally, there is a relation between classical and quantum action (transformation law) 
which reads (for x > x m i n ) 

2m(V(x) - E gr ) = 2m{V{x) - V mm ) - | ^ 2m (^)-Knm) ggn(x _ ^ (g) 

2 pm{V{x) - V min ) 

Here x m i n denotes the position of the minimum of the quantum potential. Let us 
make an ansatz for the quantum action, characterized by a mass rh and a quantum 
potential of the form 

V(x) = v 2 x 2 + v- 2 x~ 2 + v , (10) 

and see if it satisfies transformation law, Eq.([9*|). Using Eqs. (j4ll0j) . inserting it into 
Eq.([9*|). and comparing the coefficients of the terms x 2 , x~ 2 and x°, one obtains 

m v 2 = ^m 2 u 2 , rh V- 2 = ^ 2 + l} 2 , E gr = % [1 + 7] . (11) 

For the case when the classical potential is given by the parameters m — 1, h — 1, 
oo = 1 and g = 1, this yields rh v 2 — > 0.5 and rh v_ 2 — > 2. The numerical solution, 
shown in Fig. [2] gives rh v 2 —>■ 0.4999 and rh {S_ 2 — > 2.008. We observe that analytical 
and numerical results agree well. 

Dependence on boundary points. Next we look at the dependence of the quantum 
action parameters on the boundary points used in the fit. We have studied this in 




Figure 2: Inverse square potential. Classical parameters m — 1, v-i — 0.5 [to — 1), 
w_2 = 1 (g = 1). Boundary points of transition: Xi - 2 points in interval [4,5]; Xf - 
10 points in interval [0.5, 3] and At = 2 10~ 3 . 

two ways: (a) We kept the initial points Xi fixed and varied the final points Xf. (b) 
We kept the final points Xf fixed and varied the initial points X{. In the first case, 
we have considered one initial point Xi = 0.3 (kept fixed). As final points Xf we 
took 100 points uniformly distributed in an interval and varied those intervals [2,3], 
[5,6], [9, 10] and [2, 10]. In the regime T > 2 (T > 5T SC ) there is no dependence on 
the location of the final boundary points. However, for T < 2 (T < 5T SC ) there is a 
noticeable dependence. 

In the second case, as final boundary points ij we took 100 points uniformly 
distributed in the interval [2, 3] (interval kept fixed). As initial points Xi we took one 
point and varied it between 0.1 to 0.5. Those results are qualitatively the same as in 
the first case. 

One may ask: What is the reason for such dependence on boundary points? In 
our opinion the most likely explanation is that the quantum action is not an exact 
parametrisation of the transition amplitude in this regime. Hence the parametrisation 
depends on the parameters of the fit. 

The two previous cases represent a strong imbalance between the number of initial 
versus final boundary points. As a third case we have considered a more balanced 
selection of bondary points. We took 10 initial points, uniformly distributed in the 



interval [1.5,2.5] and 10 final points, uniformly distributed in the interval [1.1,2.1]. 
The global maximum relative error in fitting the transition matrix elements is about 
0.001 in the regime 1 < T < 2. 

Dependence on temporal resolution. The determination of the parameters of the 
quantum action by a global best fit requires to solve the equation of motion from S 
for all pairs of boundary points and to compute the value £ of the quantum action 
along such trajectory. Such numerical calculations depend on the temporal resolution 
At. We have studied the convergence of the parameters of the quantum action as 
a function of the density of meshpoints. We parametrize the density of meshpoints 
by taking N t meshpoints per unit time interval AT = 1. We have varied N t from 
200 up to 8000. For larger T, the results diverge. If one desires convergence for, say 
T = 8 ~ 20T SC , one needs to double the density of meshpoints. Such behavior persists 
if we want to maintain convergence for even larger T. E.g., for T — 14 w 35T SC , we 
need a 20-fold higher density of meshpoints (N t = 8000). Such rapid increase in the 
density of meshpoints could signal chaos caused by imperfect numerical solutions, 
finite internal precision, rounding errors etc. 

3. WHAT MEANS QUANTUM CHAOS? 

Many systems in nature exhibit chaos (i.e. sensitivity and exponential dependence on 
initial conditions), but their quantum analogue does not. The phenomenon of chaos 
is well understood and characterized in classical physics in terms of Poincare sections 
and Lyapunov exponents. A system is said to be chaotic if it has positive Lyapunov 
exponents or equivalently their Poincare section contains randomly distributed points. 

In general, there is no direct relation between classical and quantum dynamics ( 
Thommen et al., 2003). Some interesting behaviour is seen in systems in the semiclas- 
sical limit. For example, periodic orbits in the classical systems are very important in 
organizing the quantum eigenvalues and eigenfunctions. Contrary to this when we are 
far from the semiclassical limit the phase space notion of classical mechanics cannot 
be used in quantum physics. Heisenberg's uncertainty relation AqAp > K/2 does not 
allow to define simultanously a single point in the phase space of particle motion. 
Even though the generalized basis is free from this difficulty the sensitivity and expo- 
nential dependence on initial conditions is absent due to linearity of the Schrodinger 
equation. Necessary condition to have chaotic behavior in classical physics is the non- 
linearity introduced by dynamical variables (q,p). In quantum physics the dynamical 
variable corresponds to the wave function, whereas the potential does not depend on 
the dynamical variable. 

The quantum action is free from these difficulties. First the quantum action 
presents quantum mechanical transition amplitude by a set of parameters m and V. 



The quantum action has the same stracture as classical action, consequently there 
is no difficulty to define the phase space as in classical mechanics. Interestingly the 
quantum potential at large time limit (T — > oo) is presented by the ground state 
wave function, which means the dynamics in the quantum potential is defined by the 
dynamical variable of the Schrodinger equation. 



Let us consider as an example the following classical system (Dahlquist & Russ- 
berg, 1990), 



with the potential 

V(x,y) = v 2 (x 2 + y 2 ) + v 22 x 2 y 2 , with m = 1, v 2 = 0.5, v 22 = 0.05. (13) 

This system is globally chaotic having small islands of regularity. To find the corre- 
sponding quantum action we have parametrezed the quantum potential by analysing 
the symmetries of the original Schrodinger equation. The time-reversal symmetry, 
parity conservation and symmetry exchange x <-> y allows one to parametrize the 
quantum potential as 



where v 2 = 0.504, v 22 = 0.05 and t> 4 ~ 10~ 5 . 

Numerical studies (Caron et al.,2003) have shown the following behavior: For 
small v 22 , Poincare sections of classical and quantum systems are similar. With 
increase of energy differences between classical and quantum phase space become 
prononouced. This can be seen in Fig. where the ratio of the chaotic phase space 
volume to total phase space volume is presented for different values of energy. From 
the discussed example follows that one can define the quantum chaos as (Kroger, 
2002): Consider a classical system with action S . The corresponding quantum system 
displays quantum chaos, if the quantum action S in the large time limit T — > oo 
generates a chaotic phase space. 



We have presented the quantum action as a useful tool to study chaotic phenomenon 
in quantum systems. The structure of the quantum action allows to introduce the 
important notion of the phase space as in classical mechanics for quantum systems. 




(12) 



V(x, y) = v 2 (x 2 + y 2 ) + v 22 x 2 y 2 + w 4 (V + y 4 ) , 



(14) 



4. CONCLUDING REMARK 
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Figure 3: The estimate of the chaotic volume of the phase space. R is the ratio of the 
numeber of initial conditions with positive Lyapunov exponents to the total number 
of initial conditions. Anharmonic term t> 2 2 = 0.05 
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